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OPTICALLY  GUIDED  LASER  WAKEFIELD  ACCELERATION 

I.  INTRODUCTION 


The  possibility  of  utilizing  the  fields  of  an  intense  laser  beam  to  accelerate  particles  to 
high  energies  has  attracted  a  great  deal  of  interest  [1].  The  study  of  laser  driven  accelerators 
is  motivated  by  the  ultrahigh  fields  associated  with  high  intensity  laser  pulses.  The  peak 
amplitude  of  the  transverse  electric  field  of  the  laser  pulse  is  given  by 

£x[TV/m]  =  2.7  x  10“9/1/2[W/cm2]  =  3.2a0/A[pm],  (1) 

where  I  is  the  laser  intensity,  A  is  the  laser  wavelength,  a0  is  the  laser  strength  parameter, 
and  linear  polarization  has  been  assumed.  The  laser  strength  parameter  is  related  to  the 
power,  P,  of  a  linearly  polarized  laser  pulse  by  P[GW]  =  21.5(aoro/A)2,  where  r0  is  the 
minimum  spot  size  of  the  Gaussian  radial  profile.  Physically,  ao  is  the  normalized  trans¬ 
verse  quiver  momentum  of  the  electrons  in  the  laser  field,  a 0  =  p±/mec  =  jvj Jc,  where 
7  is  the  relativistic  factor  and  uj.  is  the  transverse  electron  velocity.  Recent  developments 
in  compact  laser  systems  [2,3]  have  resulted  in  laser  pulses  with  ultrahigh  powers,  P  ^  10 
TW,  ultrahigh  intensities,  I  >  1018  W/cm2,  moderate  pulse  energies,  W  >,  1  J,  and  ultra- 
short  pulse  durations,  77,  <  1  ps.  For  a  laser  wavelength  of  A  =  1  pm,  I  ~  1018  W/cm2 
implies  a0  ~  1  and,  hence,  highly  relativistic  electron  motion.  The  transverse  electric  field 
associated  with  this  intensity  is  El  ~  3  TV/m.  Various  laser  acceleration  concepts  involve 
transforming  a  small  fraction  of  the  ultrahigh  transverse  laser  field  into  an  axial  electric 
field  which  can  be  used  to  accelerate  electrons. 

A  laser  driven  accelerator  that  has  a  number  of  attractive  features  is  the  laser  wakefield 
accelerator  (LWFA)  [4-11].  In  the  LWFA,  a  short  (  <  1  ps)  intense  (a0  >,  1)  laser  pulse 
propagates  through  an  underdense  plasma,  A2/A2  <SC  1,  where  Ap  =  2nc/u>p  is  the  plasma 
wavelength,  up  —  {Aire2nolme)x^  is  the  plasma  frequency,  and  no  is  the  ambient  plasma 
density.  The  ponderomotive  force  associated  with  the  laser  pulse  envelope,  Fp  ~  VEl 
expels  electrons  from  the  region  of  the  laser  pulse.  If  the  laser  pulse  is  sufficiently  intense, 
virtually  all  of  the  plasma  electrons  will  be  expeiied.  When  the  laser  pulse  length,  L  =  ctl 
(defined  as  twice  the  full  width  at  half  maximum  of  the  intensity),  is  approximately  equal 
to  the  plasma  wavelength,  L  ~  Ar.  large  amplitude  plasma  waves  (wakefieids)  will  be 
excited  with  phase  velocities  approximately  equal  [12]  to  the  laser  pulse  group  velocity. 
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The  axial  and  transverse  electric  fields  associated  with  the  wakefield  can  accelerate  and 
focus  a  trailing  electron  beam.  The  ratio  of  the  accelerating  field,  Ez,  to  the  laser  field  in 
the  LWFA  is  given  by  [6,7] 

EzjEL  ~  1.2  x  10-11no'/2[cm_3]A[/im]ao(l  +  Qo/2)~1/2,  (2) 

where  linear  polarization  has  been  assumed.  For  no  =  1018  cm-3,  ao  =  1  and  A  =  1  ;tm, 
the  accelerating  field  is  Ez  ~  30  GV/m,  which  is  ~  1%  of  the  laser  field.  In  the  LWFA, 
the  plasma  serves  a  dual  purpose,  it  acts  as  a  medium  which  (i)  transforms  a  fraction  of 
the  laser  field  into  an  accelerating  field  and  (ii)  modifies  the  refractive  index  to  optically 
guide  the  laser  pulse. 

In  the  absence  of  optical  guiding,  the  interaction  distance  will  be  limited  by  diffraction 
[5].  In  vacuum,  the  laser  spot  size,  Tl ,  evolves  according  to  rL  =  ro(l  +  z2/Z^)1/2,  where  2 
is  the  axial  propagation  distance,  z  =  0  is  the  focal  point,  and  Zr  —  tttI/X  is  the  Rayleigh 
length.  Diffraction  limits  the  interaction  distance  to  ~  7 xZr.  The  maximum  energy  gain  of 
the  electron  beam  in  a  single  stage,  assuming  vacuum  diffraction,  is  AW  ~  nZREz  which, 
in  the  limit  0%  1,  may  be  written  as  AW  ~  580(A/Ap)P,  where  AW  is  in  MeV  and  P 

is  in  TW  [7].  As  an  example,  consider  a  Tl  —  200  fs  linearly  polarized  laser  pulse  with 
P  =  10  TW,  A  =  1  /urn  and  r0  =  30  pm  (ao  =  0.72).  The  requirement  that  L  ~  Ap  implies 
a  plasma  density  of  n0  =  3  x  1017  cm-3.  The  wakefield  amplitude  is  Ez  =  10  GV/m  and 
the  interaction  length  is  0.9  cm.  Hence,  a  properly  phased  tailing  electron  bunch  could 
gain  an  energy  of  AW  =  90  MeV  in  a  single  stage  without  optical  guiding. 

The  interaction  length  and  consequently  the  electron  energy  gain,  may  be  greatly 
increased  by  optically  guiding  the  laser  pulse  in  the  plasma.  Optical  guiding  in  plasmas 
can  be  achieved  by  relativistic  effects  [5-11,13-16],  preformed  density  channels  [7-9,17]  or 
tailored  laser  pulse  profiles  [8-10].  Numerical  simulations  of  three  LWFA  configurations 
which  utilize  optical  guiding  will  be  presented:  (i)  a  channel-guided  LWFA  [9],  (ii)  a 
tailored-pulse  LWFA  [10],  and  (iii)  a  self-modulated  LWFA  [11]. 

The  remainder  of  this  paper  is  organized  as  follows.  In  Section  II,  a  nonlinear,  rel¬ 
ativistic  fluid  model  is  presented  which  describes  the  intense  laser-plasma  interaction. 
This  model  is  valid  for  arbitrary  laser  intensities  and  is  capable  of  simulating,  in  the  ax- 
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isymmetric  limit,  laser  pulse  propagation  over  large  distances  (many  Rayleigh  lengths). 
Various  optical  guiding  mechanisms  are  briefly  described  in  Section  III,  including  relativis¬ 
tic  optical  guiding,  long  pulse  self-modulation,  tailored  pulse  propagation,  and  preformed 
plasma  density  channels.  In  Section  IV,  numerical  simulations  of  a  channel-guided  LWFA, 
a  tailored-pulse  LWFA,  and  a  self-mo dulated  LWFA  are  presented.  A  summary  is  given 
in  Section  V. 


II.  NONLINEAR  FORMULATION 

The  large  differences  between  the  laser  wavelength  and  other  characteristic  longi¬ 
tudinal  lengths  in  the  system,  i.e.,  laser  propagation  distance,  laser  pulse  length  and 
plasma  wavelength,  make  the  direct  numerical  integration  of  the  dynamical  equations 
over  extended  distances  impractical.  In  the  following,  a  fully  nonlinear,  relativistic,  two- 
dimensional  axisymmetric  laser-plasma  propagation  model  is  formulated  and  numerically 
evaluated  for  laser  pulses  of  ultrahigh  intensities  and  arbitrary  polarizations  [8-11].  The  for¬ 
mulation  has  a  number  of  unique  features  which  allow  numerical  simulations  to  be  carried 
out  over  extended  laser  propagation  distances.  The  appropriate  Maxwell-fluid  equations 
are  recast  into  a  convenient  form  by  (i)  performing  a  change  of  variables  to  the  “speed  of 
light”  coordinates,  (ii)  applying  the  quasi-static  approximation  (QSA)  [6],  (iii)  expanding 
in  two  small  parameters  (which  are  independent  of  the  laser  intensity)  and  (iv)  averaging 
over  the  fast  spatial  scale  length,  i.e.,  the  laser  wavelength. 

The  plasma  is  modelled  using  relativistic  cold  fluid  equations.  The  fields  associated 
with  the  intense  laser-plasma  interaction  may  be  described  by  the  normalized  scalar  and 
vector  potentials,  4>  =  e$/mec2  and  a  =  eA/raec2,  respectively.  In  the  following,  the 
Coulomb  gauge  is  used,  V  •  a  =  0,  the  ions  are  assumed  stationary  and  thermal  effects  are 
neglected.  It  is  convenient  to  introduce  the  normalized  fluid  momentum  u  =  yv/c  and  the 
fluid  quantity  p  =  n/(nory),  where  v  is  the  electron  fluid  velocity,  7  =  (1  +  u2)1/2  is  the 
relativistic  factor  of  the  electron  fluid,  and  n  is  the  electron  fluid  density.  In  the  fluid  limit, 
the  electron  current  is  given  by  J  =  — env  and  the  normalized  wave  equation  is  given  by 


1  d2\ 
c2  dt2  ) 


\d_ 
c  dt 


V0, 


(3) 
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where  kp  =  uip/c.  Poisson’s  equation  is  given  by  V20  =  k2{^ p  —  p^),  where  =  n^/n o 
and  n<°)(r)  is  the  initial  plasma  density  profile  (prior  to  the  laser  interaction)  which  can 
be  a  function  of  radial  position.  In  the  cold  fluid  limit,  the  electron  response  is  determined 
by  the  the  continuity  equation,  d(py)/dt  +  cV  •  (pu)  =  0,  and  the  momentum  equation, 
which  can  be  written  in  the  form 

1  ^  2 

-•r-(^u)  =  V(o!>  -  7)  + -ii  x  (V  x  <5u) ,  (4) 

c  at  7 

where  6u  =  u  —  a.  Equation  (4)  also  implies  that  the  quantity  h  =  V  x  6u  is  zero  provided 
it  is  initially  zero  [16].  Hence,  the  last  term  on  the  right  side  of  Eq.  (4)  will  be  neglected. 

The  Maxwell-fluid  equations  can  be  simplified  by  the  following.  The  full  set  of  equa¬ 
tions  is  recast  into  speed  of  light  coordinates  by  introducing  the  independent  variables 
C  =  z  —  ct  and  r  =  t.  The  QSA  is  then  applied  [6].  In  the  QSA,  the  electron  transit  time 
through  the  laser  pulse  (~  tl)  is  assumed  to  be  short  compared  to  the  laser  pulse  evolution 
time,  which  is  determined  by  the  pulse  diffraction  time,  ZR/c,  or  by  the  pulse  dispersion 
time,  u)/up,  where  ui  is  the  laser  frequency.  In  the  QSA,  the  electrons  experience  essentially 
static  fields,  allowing  the  djdr  derivatives  to  be  neglected  in  the  fluid  equations,  but  not 
in  the  wave  equation.  Within  the  QSA,  the  quantity  7  —  uz  —  <p+az  is  an  invariant  which  is 
set  equal  to  unity,  i.e.,  its  value  prior  to  the  arrival  of  the  laser  pulse.  Hence,  7— uz  —  1  +  0, 
where  0  =  <f>  —  az  is  the  wake  potential.  The  resulting  equations  are  expanded  to  first 
order  in  the  parameters  ei  =  1/kri  «C  1  and  62  =  kp/k  1,  where  k  —  2tt/X.  All  the 
fluid  and  field  quantities  are  expanded  in  slow  and  fast  terms,  i.e.,  Q  =  Q„  +  Qf,  where 
\dQf/d£\  ~  |fcQ/|  and  \dQ„/d£\  ~  jfcpQ,j.  Within  this  representation,  the  nonlinear  fluid 
equations  are  averaged  over  the  laser  wavelength  in  the  C, r  frame.  The  ^-averaging  allows 
for  all  the  laser-plasma  response  quantities  to  be  evaluated  on  the  slow  spatial  scale,  i.e., 
A p  or  L,  permitting  solutions  over  extended  propagation  distances. 

The  resulting  equations  describe  the  slowly  varying  components  of  the  fluid  and  field 
quantities  [8-11]: 

V^a  =  k^pu  -  (5a) 

(vl  +  |^0  =  ^(7  P-P(0)),  (56) 


4 


~  (u  -  a)  =  V  (7  -  0) ,  (5c) 

where  the  subscript  s,  denoting  the  slow  component  of  the  quantity,  has  been  dropped. 
Equations  (5a-c)  represent,  respectively,  the  slow  components  of  the  wave,  Poisson’s,  and 
momentum  equations.  The  slow  component  of  the  continuity  equation  follows  from  the 
divergence  of  Eq.  (5a).  The  slow  component  of  the  relativistic  factor  is 

_  l  +  uj  +  |Q/l2/2  +  (l  +  ^)2 

2(1  + VO  ‘  { 

In  Eq.  (6)  and  throughout  the  remainder  of  this  paper,  a  linearly  polarized  laser  pulse 
with  amplitude  |d/|  is  assumed.  The  transverse  component  of  the  laser  radiation  field 
is  a/  =  a/(r,  £,  r)  exp(i&C)/2+  c.c.,  where  a f  is  the  complex,  slowly  varying  amplitude, 
which  satisfies  the  wave  equation 


Typically,  the  last  term  on  the  left  of  Eq.  (7)  is  small  and  can  be  neglected,  as  is  done 
in  the  simulations  discussed  below.  Within  the  QSA,  the  self-consistent,  slowly  varying 
equations  in  the  £,  r  variables,  describing  the  laser-plasma  interaction,  to  first  order  in  ei 
and  C2i  are  given  by  Eqs.  (5)-(7). 

Equations  (5a)-(5c)  can  be  combined  to  yield  a  single  equation  for  ip  in  terms  of  jd/|2  of 
the  form  d2ip/d(,2  —  G(ip,  |d/|2),  where  G  is  an  involved  function  [8-11].  The  equation  for 
ip  is  obtained  by  noting  that  (1  +  ip)p  =  p +  k~2V\ip  and  k2pu±  =  V±dip/d(p,  assuming 
axisymmetry.  Note  also  that  the  refractive  index  is  solely  a  function  of  ip  through  p,  i.e., 
r]R  =  1  —  k2p/2k 2.  Equation  (7)  together  with  d2ipjdC,2  =  G  completely  describe  the 
2D-axisymmetric,  quasi-static  laser-plasma  interaction.  In  this  axisymmetric  model,  the 
effects  of  various  instabilities,  such  Raman  sidescatter  and  stimulated  backscatter,  will  be 
neglected  [15,18,19].  The  wake  potential,  ip,  is  related  to  the  axial  electric  field,  Ez,  of 
plasma  response  (wakefield)  by  Ez  =  -dip/d(,  where  Ez  —  eEz/mec2.  Furthermore,  Eqs. 
(5)-(7),  in  the  appropriate  limits,  reduce  to  models  which  have  been  previously  studied.  In 
the  broad  pulse  limit,  V±  —> >  0,  Eqs.  (5)-(7)  reduce  to  a  set  of  coupled  equations  describing 
the  well-know  one-dimension  model  [6,20].  In  the  limit  djdC,  — »  9,  Eqs.  (5)-(7)  reduce  to 
a  single  equation  describing  the  evolution  of  a  long,  axially  uniform  pulse  [14]. 
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III.  OPTICAL  GUIDING  MECHANISMS 


The  optical  guiding  mechanisms  discussed  below  are  based  on  the  principle  of  refrac¬ 
tive  guiding.  Refractive  guiding  becomes  possible  when  the  radial  profile  of  the  index  of 
refraction,  r?R(r),  exhibits  a  maximum  on  axis,  i.e.,  df]Rjdr  <  0.  The  general  expression 
for  the  index  of  refraction  for  a  large  amplitude  electromagnetic  wave  in  a  plasma  is  given 
by  *?H  —  1 —  wpP(r)/2w2,  where  p  =  n(r)/n07(r).  Hence,  the  r)R(r )  profile  can  be  modified 
by  (i)  the  relativistic  transverse  quiver  motion  of  the  electrons  in  the  laser  field,  7  =  7±(r), 
where  7x  =  [1  +  |a./(r)|2/2]1/2T  (ii)  a  preformed  plasma  density  channel,  n  =  n(°)(r),  and 
(iii)  the  effects  of  a  plasma  wave,  n  —  no  -f  8n(r),  where  6n  —  6no(r)  exp  ikp£  is  the  den¬ 
sity  oscillation.  These  three  effects  are  responsible  for  relativistic  optical  guiding,  density 
channel  guiding,  and  the  self-modulation  of  long  laser  pulses. 

A.  Relativistic  Optical  Guiding 

In  the  standard  theory  of  relativistic  optical  guiding  [5,13,14],  only  the  effects  of 
the  transverse  quiver  motion  of  the  electrons  are  included  in  the  expression  for  r]R,  i.e., 
n  =  no  and  7  —  y±(r).  Analysis  of  the  wave  equation  with  an  index  of  refraction  of 
this  form  indicates  that  when  the  laser  pulse  power  exceeds  a  critical  power,  P  >  P criti 
where  Peru  —  17{(jj/up)2  GW,  relativistic  effects  can  prevent  the  diffraction  of  the  laser 
pulse.  The  self-consistent,  nonlinear  theory  developed  in  Refs.  [6,9]  showed,  however,  that 
relativistic  optical  guiding  is  ineffective  in  preventing  the  diffraction  of  sufficiently  short 
pulses,  L  <  Ap/(1  +  |d/|2/2)1//2.  This  is  due  to  the  fact  that  the  index  of  refraction  becomes 
modified  by  the  laser  pulse  on  the  plasma  frequency  time  scale,  not  the  laser  frequency 
time  scale.  Typically,  relativistic  guiding  only  effects  the  body  of  long  pulses,  L  >  Xp. 

B.  Long  Pulse  Self-Modulation 

Recently,  it  was  shown  that  long  pulses,  L  >  Ap,  with  P  >  Pcrit  can  be  subject  to 
severe  self-modulation  [8,9,11].  Specifically,  a  plasma  wave,  excited  by  either  the  finite  rise 
of  the  laser  pulse  envelope  or  by  a  forward  Raman  scattering  (FRS)  instability,  can  strongly 
affect  the  focusing  properties  of  the  pulse  body.  Consider  a  long  laser  pulse,  L  3>  Ap,  with 
P  =  Pcrit  on  which  a  finite  wakefield  (a  plasma  wave  with  phase  velocity  approximately 
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equal  to  the  group  velocity  of  the  laser  pulse)  exists.  The  wakefieid,  which  consists  of 
a  plasma  density  modulation  of  the  form  6n  =  6n0(r)  cos(kp(),  modifies  the  plasma’s 
refractive  index  [17].  In  regions  of  a  local  density  channel,  i.e.,  where  d bnjdr  >  0,  the 
radiation  focuses.  In  regions  where  ddn/dr  <  0,  diffraction  is  enhanced.  This  causes  the 
laser  pulse  envelope  to  become  modulated  at  the  plasma  wavelength,  which  subsequently 
enhances  the  growth  of  the  plasma  wave,  and  the  process  proceeds  in  a  highly  nonlinear 
manner.  In  the  limit  r*o  >  Ap  and  P  <  Pent,  in  which  the  e-folding  length  of  the  standard 
FRS  instability,  Lrrsi  is  much  less  than  the  Rayleigh  length,  i.e.,  Lrrs  <C  Zr,  the  self¬ 
modulation  process  is  adequately  described  by  standard  FRS  theory.  In  the  limit  P  >  Pent 
and  Lrrs  £  Zr,  however,  the  modulation  process  is  highly  nonlinear  with  a  growth  rate 
much  greater  than  that  expected  from  standard  FRS  theory  [11].  This  results  in  a  fully 
self-modulated  laser  pulse,  composed  of  a  series  of  laser  “beamlets”  of  length  ~  Ap/2, 
which  can  remain  optically  guided  over  several  Rayleigh  lengths.  Associated  with  the 
periodic  beamlet  structure  axe  large  amplitude  wakefields  which  can  trap  and  accelerate  a 
trailing  electron  beam.  This  process  forms  the  basis  of  the  self-modulated  LWFA  discussed 
below  [11]. 

C.  Tailored  Pulse  Propagation 

A  laser  pulse  with  an  appropriately  tailored  envelope  can  propagate  many  Rayleigh 
lengths  without  significantly  altering  its  original  profile  [8-10],  Consider  a  long  laser  pulse, 
L  »  Ap,  in  which  the  spot  size  is  tapered  from  a  large  value  at  the  front  to  a  small  value  at 
the  back,  so  that  the  laser  power,  P  ~  r2  )d/}2,  is  constant  throughout  the  pulse  and  equal 
to  Pcrit  •  The  leading  portion  (<§:  Ap)  of  the  pulse  will  diffract  as  if  in  vacuum  [6,8].  However, 
since  r*,  is  large  at  the  front  of  the  pulse,  the  Rayleigh  length  is  also  large.  Hence,  the 
locally  large  spot  size  allows  the  pulse  front  to  propagate  a  long  distance,  whereas  the  body 
of  the  pulse  will  be  relativistically  guided.  Also,  since  |a/|2  increases  slowly  throughout 
the  pulse,  detrimental  wakefieid  effects  (i.e.,  self-modulation)  are  reduced  [8-10,21]. 

D.  Plasma  Density  Channel  Guiding 

A  preformed  plasma  density  channel  can  guide  short,  intense  laser  pulses  [7-9,17], 
In  the  weak  laser  pulse  limit,  |d/|2  1,  the  index  of  refraction  is  given  by  r?R  =  1  - 


7 


ujpP^ /2uj2 .  Optical  guiding  requires  drfRjdr  <  0,  hence,  a  preformed  density  channel, 
n(°)(r)  =  noP(°Hr)i  can  prevent  pulse  diffraction.  Analysis  of  the  wave  equation  in  the 
weak  pulse  limit  indicates  that  a  parabolic  density  channel  will  guide  a  Gaussian  laser 
beam  provided  that  the  depth  of  the  density  channel  [7-9,17]  is  An  =  l/7rrer£,  where 
An  =  u/'°)(r£,)  —  n^(0)  and  re  -  e2 jmec2  is  the  classical  electron  radius.  One  possible 
method  for  creating  a  preformed  density  channel  is  to  use  the  ponderomotive  force  of  one 
or  more  additional  laser  beams. 

IV.  LWFA  SIMULATIONS 

To  study  laser  pulse  propagation  and  wakefield  generation  in  plasma,  the  set  of 
Maxwell  fluid  equations,  Eq.  (5)-(7),  are  solved  numerically.  Specn'.eally,  the  equation 
d2%j;/d£2  =  G(^’,|d/|2)  is  solved  for  a  given  l»~er  envelope,  and  the  wavi  .quation,  Eq. 
(7),  is  solved  to  update  the  laser  envelope,  neglecting  the  d2fdr 2  term.  Retention  of  the 
d2 /d^dr  term  Eq.  (7)  is  necessary  to  correctly  describe  the  group  velocity,  vg  <  c,  of  the 
laser  pulse,  the  phase  velocity  of  the  wakefield,  and  the  phase  detuning  of  the  acceler¬ 
ated  electron  beam  [11,12].  To  study  the  acceleration  and  trapping  of  electrons  by  the 
high-gradient  wakefield,  a  particle  code  is  used  to  evolve  a  distribution  of  non-interacting 
test  electron  beam  particles  in  the  time-resolved  electric  and  magnetic  wakefields  of  the 
simulation.  Examples  of  three  LWFA  acceleration  schemes  will  be  presented,  a  channel- 
guided  LWFA  [9],  a  tailored-pulse  LWFA  [10],  and  a  self-modulated  LWFA  [11].  In  the 
channel-guided  LWFA  and  the  tailored-  pulse  LWFA  simulations,  in  which  phase  detuning 
is  not  an  issue,  the  term  proportional  to  d2 /d^dr  has  also  been  neglected  on  the  left  side 
of  the  wave  equation,  Eq.  (7). 

A.  Channel-Guided  LWFA 

The  results  of  a  channel-guided  LWFA  simulation  [9]  are  shown  in  Figs.  1-3.  In  this 
example,  the  initial  axial  laser  profile  is  given  by  |d/(0!  =  a0sm(-n^/L)  for  0  <—(<£, 
with  a0  =  0.72  and  L  =  120  pm  (400  fs).  Also,  A  =  1  pm  and  rL  =  60  ^m  (Gaussian 
radial  profile),  which  implies  Zr  =  1.1  cm  and  P  -  40  TW.  The  density  on  axis  is  chosen 
such  that  L  =  Xp  (no  =  7.8  x  1016  cm-3)  and  a  preformed  density  channel  with  a  parabolic 
profile  is  assumed  with  An  =  l/7rrer£  =  3.2  x  1016  cm-3. 
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Figure  1(a)  shows  the  evolution  of  the  laser  spot  size  versus  propagation  distance,  cr. 
The  laser  pulse  remains  well-guided  by  the  density  channel,  the  laser  spot  size  exhibiting 
small  oscillations  about  its  initial  value  over  the  full  20 Zr  —  23  cm  simulation  length. 
After  cr  =  20 Zr,  the  pulse  profile  shows  very  little  distortion  from  its  initial  profile.  A 
surface  plot  of  the  electron  density  profile  at  cr  =  20 Zr  is  shown  in  Fig.  2.  The  initial, 
unperturbed  parabolic  profile  can  be  seen  at  the  right  (£  =  0),  and  the  distortion  of  the 
channel  by  the  laser  pulse,  including  the  excitation  of  a  large  amplitude  wakefield  along 
the  axis,  is  evident  in  the  region  £  <  0.  In  this  example,  nearly  all  the  electrons  have  been 
expelled  from  the  vicinity  of  the  laser  pulse.  A  plot  of  the  axial  electric  field  along  the  axis 
at  cr  =  20 Zr  is  shown  in  Fig.  3.  The  effects  of  the  wakefields  on  a  continuous  electron 
beam  with  an  initial  normalized  emittance  e„  =  1.0  mm-mrad,  RMS  radius  r*>  =  10  p m 
and  energy  Et,  =  2  MeV  was  also  simulated.  After  cr  =  20  cm,  approximately  70%  of  the 
beam  electrons  were  trapped  and  accelerated.  The  peak  electron  energy  increases  nearly 
linearly  with  propagation  distance  with  an  average  acceleration  gradient  of  5.25  GeV/m 
(1  GeV  in  20  cm). 


B.  Tailored-Pulse  LWFA 

Simulation  results  for  a  LWFA  utilizing  a  tailored  laser  pulse  [10]  are  shown  in  Figs.  1,4 
and  5.  The  initial  normalized  intensity  profile,  |d/|2,  is  shown  in  Fig.  4(a).  The  local  spot 
size  at  the  front  (£  =  0)  of  the  pulse  is  large,  tr  ~  8AP,  and  tapers  down  over  £  =  — 2AP  to 
rL  =  Ap  (a  Gaussian  radial  profile  assumed  throughout).  The  initial  axial  laser  envelope  is 
given  by  [&/(£)]  =  sin(— 7r£/4Ap)  for  0  <  — £  <  2Ap  such  that  P  =  Pcrit  at  each  £  slice, 
i.e.,  r£,(£)ao(£)  =  0.9Ap.  Also,  A  =  1  p m  and  Ap  =  30  pm  (no  =  1.2  x  1018  cm-3,  initially 
a  uniform  plasma),  such  that  P  —  PCTit  —  16  TW.  This  gives  a  peak  value  of  a0  =  0.9  at 
the  back  of  the  pulse  where  ri  =  Ap,  which  corresponds  to  a  Rayleigh  length  of  Zr  =  0.28 
cm.  The  pulse  intensity  then  terminates  over  a  distance  of  Lfau  =  Ap/2.  The  pulse  energy 
is  approximately  3  J  and  the  pulse  length  is  approximately  L  =  2AP  =  60  pm  (200  fs). 
Because  Lfaii  <  Ap,  a  large  amplitude  wakefield  will  be  excited  behind  the  pulse. 

Figure  4(b)  shows  the  normalized  intensity  profile  after  propagating  cr  =  10 Zr  =  2.8 
cm.  The  pulse  is  somewhat  distorted,  but  largely  intact.  The  evolution  of  the  pulse  spo* 
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size  at  the  position  of  peak  intensity  verses  propagation  distance  is  shown  in  Fig  1(b), 
indicating  that  guiding  has  been  achieved  over  the  10 Zr  simulation  region.  The  axial 
electric  field  of  the  wake  on  axis  after  cr  =  10 Zr  is  shown  in  Fig.  5.  The  evolution 
of  a  continuous  electron  beam  with  an  initial  normalized  emittance  £„  =  10  mm-mrad, 
RMS  radius  rb  =  5  /an  and  energy  Eb  =  2  MeV  was  simulated  using  the  self-consistent 
wakefields.  After  cr  =  2.8  cm,  approximately  60%  of  the  beam  electrons  were  trapped  and 
accelerated.  The  peak  energy'  of  tue  beam  electrons  experienced  an  average  gradient  of  27 
GeV/m  (750  MeV  in  2.8  cm). 


C.  Self-Modulated  LWFA 

The  mechanism  of  long-pulse  self-modulation  can  form  the  basis  for  a  LWFA  Ill). 
To  operate  in  the  self-modulated  regime,  it  is  necessary  that  (i)  the  pulse  length  be  long 
compared  to  the  plasma  wavelength,  L  >  Ap,  and  (ii)  the  power  be  equal  or  greater  than 
the  critical  power,  P  >  PCTit ■  For  fixed  laser  pulse  parameters,  i.e.,  fixed  L  and  P,  the  self- 
modulated  regime  can  be  reach  by  operating  at  sufficiently  high  density,  since  Ap  ~  1  /y/no 
and  PCrit  ~  l/n0-  The  advantages  of  the  self-modulated  LWFA  over  the  standard  LWFA 
(in  which  L  —  Ap)  are  simplicity  and  enhanced  acceleration.  Simplicity  in  that  a  matching 
condition  of  L  ~  Ap,  a  preformed  density  channel  and/or  special  pulse  tailoring  are  not 
required.  Enhanced  acceleration  is  achieved  for  several  reasons:  (i)  The  self-modulated 
LWFA  operates  at  higher  density,  hence  a  larger  wakefield  will  be  generated,  since  E .  ~ 
l/v'nj,  as  indicated  by  Eq.  (2);  (ii)  since  P  >  Pent,  the  laser  pulse  will  tend  to  focus 
to  a  higher  intensity,  thus  increasing  a0  and  Ez\  (iii)  the  wakefield  is  resonantly  excited, 
i.e.,  excited  by  a  series  of  beamlets  as  opposed  to  a  single  pulse  as  in  the  standard  LWFA; 
and  (iv)  relativistic  optical  guiding  allows  the  modulated  pulse  structure  to  propagate  for 
several  Rayleigh  lengths,  thus  extending  the  acceleration  distance.  The  disadvantages  of 
the  self-modulated  LWFA  are  (i)  at  higher  densities  the  laser  pulse  group  velocity  (~  the 
plasma  wakefield  phase  velocity  [12])  decreases  and,  hence,  electron  dephasing  from  the 
plasma  wakefield  can  limit  the  acceleration  distance,  and  (ii)  the  modulated  pulse  structure 
eventually  diffracts. 

The  properties  of  the  self-modulated  LWFA  are  illustrated  by  the  following  simula- 
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tions.  For  fixed  laser  pulse  parameters,  two  cases  will  be  considered:  (I)  a  standard  LWFA 
in  which  L  ~  Ap  and  P  <  Peru  and  (II)  a  self-modulated  LWFA,  in  which  L  >  Xp  and 
P  >  Pcrit •  The  laser  parameters  for  both  these  cases  are  identical:  a  Gaussian  axial  inten¬ 
sity  profile  with  a  pulse  length  L  —  90  (300  fs),  A  =  1  /xm,  ao  =  0.7,  r0  =  31  jxm  (in 

vacuum),  which  corresponds  to  Zr  =  3  mm,  P  =  10  TW  and  a  pulse  energy  of  W  =  1.5 
J,  well  within  the  bounds  of  present  technology.  The  simulation  begins  at  r  =  0  as  the 
laser  pulse  enters  the  plasma,  initially  converging  such  that  in  vacuum  it  would  focus  to  a 
minimum  spot  size  of  ro  =  31  fxm  at  cr  =  3 Zr.  The  plasma  density  is  initially  increasing, 
reaching  full  density  at  cr  =  2 Zr.  The  simulation  continues  until  cr  =  10 Zr  =  3  cm. 
In  both  cases,  the  acceleration  and  trapping  of  a  continuous  electron  beam  with  initial 
energy  of  3  MeV  and  normalized  emittance  en  =  130  mm-mrad  is  considered.  The  elec¬ 
tron  beam  is  initially  converging  such  that  in  vacuum  it  would  focus  to  a  minimum  RMS 
radius  rj,  =  200  /xm  at  cr  =  3 Zr.  With  such  a  large  initial  emittance,  only  a  small  fraction 
(~  1%)  of  the  particles  will  be  trapped  and  accelerated. 

For  the  standard  LWFA,  Case  I,  the  requirement  L  =  Ap  =  90  /xm  implies  a  density 
of  n0  =  1.4  x  1017  cm-3.  At  this  density,  P  <  Peru  =  140  TW,  such  that  relativistic 
guiding  effects  are  unimportant.  In  fact,  the  presence  of  the  plasma  has  little  effect  on  the 
evolution  of  the  laser  pulse,  which  reaches  a  peak  intensity  of  |d/j2  =  0.56  at  cr  =  3 Zr. 
The  evolution  of  the  spot  size,  Fig.  1(c),  is  very  close  to  vacuum  diffraction.  This  is  also 
evident  in  Fig.  6(a)  (dashed  line),  where  the  peak  accelerating  field,  plotted  versus  time,  is 
symmetric  about  the  focus,  cr  =  3 Zr.  After  cr  =  10 Zr  =  3  cm,  a  small  fraction  (~  0.1%) 
of  the  test  electron  beam  particles  have  been  trapped  and  accelerated.  At  cr  =  2  cm,  the 
peak  particle  energy  is  48  MeV,  which  implies  an  average  acceleration  of  2.4  GeV/m,  as 
shown  in  Fig.  7(a)  (dashed  line). 

For  the  self-modulated  LWFA,  Case  II,  the  density  is  increased  such  that  P  =  1.5Pcr,t 
=  10  TW,  which  implies  no  =  2.8  x  1018  cm-3,  which  is  20  times  higher  than  in  Case  I. 
At  this  density  L  >  Ap  =  20  /'m,  i.e.,  the  laser  pulse  now  extends  over  ~  4.5Ap.  Figure 
8  shows  the  laser  intensity  at  (a)  cr  =  2 Zr  and  (b)  cr  =  3.2 Zr.  The  axial  electric 
field  on  axis  at  cr  =  3.2 Zr  is  shown  in  Fig.  9.  The  laser  pulse  has  become  modulated 
(three  peaks  are  observable,  separated  by  Ap)  and  the  plasma  wave  is  highly  nonlinear. 
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In  addition,  relativistic  optical  guiding  effects  have  focused  the  laser  to  a  much  higher 
intensity  than  was  observed  in  Case  I.  The  evolution  of  the  laser  spot  size  is  shown  in  Fig. 
1(d),  indicating  that  the  pulse  has  focused  to  a  smaller  spot  size  and  remains  guided  over 
~  5.5Zr.  A  plot  of  the  peak  accelerating  field  versus  time,  Fig.  6(b)  (solid  fine),  shows 
that  the  highly  nonlinear  fields  persist  as  the  laser  pulse  is  optically  guided.  A  maximum 
accelerating  field  of  ~  130  GV/m  was  obtained.  Approximately  2%  of  the  test  electron 
beam  particles  were  trapped  and  accelerated.  The  peak  particle  energy  of  430  MeV  is 
observed  at  cr  =  6 Zr  =  1.8  cm.  At  cr  =  lOZ^  —  3  cm,  however,  the  peak  particle  energy 
has  dropped  to  290  MeV  due  to  the  reduced  group  velocity  of  the  laser  pulse,  which  causes 
the  electrons  to  slip  out  of  phase  with  the  wakefield  and  become  decelerated.  Figure  7(b) 
(solid  line)  shows  acceleration  to  430  MeV  over  1.8  cm  which  gives  an  average  gradient  of 
24  GeV/m.  This  is  an  order  of  magnitude  increase  compared  to  the  standard  LWFA  of 
Case  I. 


V.  CONCLUSION 

The  laser  wakefield  accelerator  [4-11]  is  capable  of  generating  ultrahigh  accelerating 
and  focusing  fields  over  extended  distances  (many  Rayleigh  lengths),  which  can  trap  and 
accelerate  a  trailing  electron  beam.  The  large  amplitude  wakefields  are  generated  by  the 
ponderomotive  forces  associated  with  the  envelope  of  the  laser  pulse.  These  forces  can 
expel  electrons  from  the  vicinity  of  the  laser  pulse,  thus  exciting  a  large  amplitude  plasma 
wave  (wakefield)  with  phase  velocity  near  the  speed  of  light.  For  a  single  laser  pulse,  the 
maximum  wakefield  amplitude  results  when  the  pulse  length  is  approximately  equal  to  the 
plasma  wavelength,  L  ~  Ap.  Diffraction  places  the  most  severe  limitation  on  the  pulse 
propagation  distance,  hence  optical  guiding  is  necessary  to  accelerate  electrons  to  high 
energies  [5].  Laser  pulse  diffraction  can  be  prevented  by  the  effects  of  relativistic  optical 
guiding,  preformed  density  channels,  and/or  appropriately  tailored  pulse  profiles  [5-11,13- 
17J.  To  study  pulse  propagation  and  wakefield  generation,  an  axisymmetric  Maxwell-fluid 
model  was  developed  [8-11].  This  model  is  valid  for  arbitrary  laser  intensities  and  can 
simulate  laser  propagation  in  plasma  over  large  distances  (many  Rayleigh  lengths). 

Three  possible  LWFA  configurations  were  simulated:  (i)  a  channel-guided  LWFA  [9], 
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(ii)  a  tailored-pulse  LWFA  [10],  and  (iii)  a  self-modulated  LWFA  [11],  The  channel- 
guided  LWFA  and  tailored-pulse  LWFA  simulations  demonstrate  that  a  laser  pulse  can  be 
optically  guided  over  many  Rayleigh  lengths  while  generating  large  amplitude  wakefields, 
which  can  trap  and  accelerate  a  test  electron  beam  to  high  energies  (~  1  GeV).  The  self- 
modulated  LWFA,  a  recently  proposed  configuration,  is  of  particular  interest,  especially 
concerning  the  design  of  proof-of-principle  LWFA  experiments  [ll].  In  the  self-modulated 
LWFA,  the  plasma  density  is  chosen  to  be  sufficiently  high  so  that  the  laser  pulse  is  long 
compared  to  a  plasma  wavelength,  L  >  Ap,  and  the  power  exceed  the  critical  power,  P  > 
Per  it-  Previously  it  was  believed  that  in  the  long-pulse  regime,  in  which  \daf  /dQ  ~  L|a/| 
and  L  Ap,  the  plasma  electrons  respond  adiabatically,  and  no  significant  wakefield  would 
be  generated.  This,  however,  is  not  necessarily  the  case.  When  L  >  Ap  and  P  >  Pcr«t> 
the  laser  pulse  can  become  strongly  self-modulated  at  the  plasma  wavelength  and  large 
amplitude  wakefields  are  resonantly  excited.  This  results  in  enhanced  wakefield  amplitudes 
and  enhanced  acceleration  in  comparison  to  the  standard  LWFA  (in  which  L  ~  Ap).  For 
the  simulation  presented  in  Section  IV,  electron  acceleration  was  enhanced  by  an  order  of 
magnitude  by  operating  in  the  self-modulated  regime  for  fixed  laser  parameters,  simply 
by  increasing  the  ambient  density  by  a  factor  of  20.  The  self-modulated  LWFA  has  the 
additional  advantage  of  simplicity,  i.e.,  preformed  density  channels,  pulse  profile  tailoring, 
and  a  matching  condition  of  L  ~  Ap  are  not  necessary.  Hence,  the  self-modulated  regime 
is  well-suited  for  a  proof-of-principle  LWFA  experiment. 
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Fig.  1.  Laser  spot  size  normalized  to  the  plasma  wavelength,  r^f  Xp,  at  the  position  of 
peak  intensity,  versus  propagation  distance  normalized  to  the  Rayleigh  length, 
ct/Zr,  for  (a)  the  channel-guided  LWFA,  (b)  the  tailored-pulse  LWFA,  (c)  vac¬ 
uum  diffraction  for  the  parameters  of  the  self-modulated  LWFA,  and  (d)  the 
self-modulated  LWFA. 
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Fig.  3.  Axial  electric  field  Ez  on  axis  at  cr  =  20 ZR  for  the  channel-guided  LWFA. 
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Fig.  7.  Peak  energy  of  the  trapped  electrons  verses  propagation  distance,  , 
standard  LWFA  (Case  I),  dashed  line;  and  (b)  the  self -modulated 
II),  solid  line. 
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